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Abstract 

We study the quantization of Hitchin systems in terms of /3-deformations of generalized matrix 
models related to conformal blocks of Liouville theory on punctured Riemann surfaces. We show 
that in a suitable limit, corresponding to the Nekrasov-Shatashvili one, the loop equations of 
the matrix model reproduce the Hamiltonians of the quantum Hitchin system on the sphere 
and the torus with marked points. The eigenvalues of these Hamiltonians are shown to be the 
ei-deformation of the chiral observables of the corresponding M = 2 four dimensional gauge 
theory. Moreover, we find the exact wave-functions in terms of the matrix model representation 
of the conformal blocks with degenerate field insertions. 



1 Introduction 



M = 2 supersymmetric gauge theories in four dimensions display very interesting mathemat- 
ical structures in their supersymmetrically saturated sectors. These structures allow an exact 
characterization of several important physical aspects, such as their low energy behavior and 
stable spectra. These data are encoded in the celebrated Seiberg-Witten solution |T]. It was 
soon realized that the Seiberg-Witten data can be recovered from integrable systems in terms 
of their spectral curves [2]. In this context the Hitchin integrable system has emerged as 
the fundamental geometric structure underlying the M-theory description of M = 2 theories 

mm 

On the other hand the Seiberg-Witten solution can also be recovered, at least in the case of 
linear and elliptic quiver M = 2 theories, via equivariant localization on the instanton moduli 
space [6j [T] . Indeed, this approach contains further information encoded in the expansion in the 
equivariant parameters of the f2-background. This opens the issue of relating the full Nekrasov 
partition function to a suitable quantization of the Hitchin system. 

A crucial result in this context is provided by the AGT correspondence [8] relating the 
Nekrasov partition function to conformal blocks of Liouville/Toda field theories in two dimen- 
sions. In [91 dm it was proposed that this correspondence should be regarded as a two parameter 
quantization of the Hitchin system itself, or, in field theory language, as its second quantization. 
Here we will address these issues in the particular limiting case in which one of the two equiv- 
ariant parameters is vanishing. This was identified by Nekrasov and Shatashvili [11] to provide 
the first quantization of the integrable system*]- In the context of AGT correspondence the 
instanton partition function in the Nekrasov-Shatashvili limit can be related to the insertion of 
degenerate fields in the Liouville theory [28l EH], which corresponds to the insertion of surface 
operators in the gauge theory side [301 EH E21 E31 EH ESI ES] . 

In our approach, we will make use of the matrix model perspective on AGT correspondence 
developed in [37]. This was further elaborated for Liouville theory on the sphere in [3H1 ES, 
SQl SH S3 SHI S3l SHI SSI SZl SHI SHI EQl EH [521 [53], on the torus in [Ml [55] and in [56] at 
all genera. In this context the equivariant parameters of the Nekrasov partition function are 
encoded in the /3-deformation of the standard Van- der- Monde measure [37]. Notice that for 
/3-deformed matrix models the algebraic equation defining the spectral curve gets deformed into 
a differential equation which can be interpreted as a Schrodinger equation [57J, EE1 EH]. Our 
proposal identifies this differential equation, in the Nekrasov-Shatashvili limit, as providing the 
quantum Hamiltonians of the associated Hitchin integrable system. Moreover the associated 



*See also [HI [131 IHl H3 HEl H3 HH1 US HOI [2H HSl |Z3 for the relation between the gauge 

theory and the quantized integrable system. 
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wave- function is described in terms of the /3-deformed generalized matrix model corresponding 
to degenerate field insertions in the Liouville theory [301 ED EH] . 

Let us notice also that the quantization of Hitchin systems plays a vital role in the context 
of Langlands duality [60] . 

The organization of this paper is as follows. In section El we review the M-theory perspective 
on the Hitchin system and some basic facts on its quantization at low genera, namely the sphere 
and the torus with marked points. We derive the loop equation for the generalized /9-deformed 
matrix model both in the sphere and torus case in sections [3]and|4]respectively and we show that 
in the Nekrasov-Shatashvili limit these reproduces the Hamiltonians of the quantized Hitchin 
system. Moreover, we provide a description of the associated wave-functions in terms of the 0- 
deformed generalized matrix model describing degenerate field insertions in the Liouville theory 
[30] |3T| [29] . In section [5] we present some final comments and further directions. 

2 Quantum Hitchin systems in nuce 

In the M-theory framework, the Hitchin system arises by considering the geometry of a system 
of N M5-branes wrapped on a manifold Yq with topology C x R 4 x {pt} in T*C x R 4 x R 3 where 
C is a Riemann surface. This [37] should be equipped with a non trivial fibration of R 4 over C 
which specifies the ^-background of Nekrasov [6]. The geometry of the M5-branes bound state 
is described by an iV-fold branched covering of C given by the algebraic equation 



i=2 

where x is a section of T*C and <pj are (j, 0)-holomorphic differentials on C whose singularity 
structure at the punctures identifies the matter content of the gauge theory. The quantum 
mechanics of this latter structure passes by interpreting (12.11) as the spectral curve of the 
associated classical system encoding the Seiberg- Witten solution (I2.ip , and then quantizing via 
a suitable deformation. 

In what follows, we focus on the case with two M5-branes and the corresponding Hitchin 
systems. In this case, the M-theory curve is specified by the quadratic differential: x 2 = ^2(2). 
We further focus on the situation where the singularities of the quadratic differential are all of 
regular type meaning the poles of at most degree 2. There is one kind of regular punctures in 
this two M5-branes case, so the Riemann surface C is just specified by genus and the number of 
punctures. Thus, we denote this by C 9i „. Under the particular marking of C 5i „, the worldvolume 
low energy theory is the weakly coupled SU(2) n+39 ~ 3 superconformal quiver gauge theory [611 H]. 
The Hitchin system is associated with this gauge theory. 



N 




(2.1) 
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In order to have explicit parametrization of the Hamiltonians, let's review and specify 
Hitchin systems and their quantization at low genera, that is sphere and torus with an arbitrary 
number of punctures, following the approach of [62] . 

Let us start from some general features of the system. The Higgs field satisfies the condition 



d A $ = h$« 

k 

which can be simplified by A = h~ 1 dh and $ = h _1 ^h to 

<9<£ = 5>*<W ( 2 - 2 ) 

k 

where Vk = h^iwkj^khiwk)- Eq. fl2.2p admits a unique solution iff J2k Uk = 0> given by 

$ = ^i/ fc a W * + $° (2.3) 

k 

where upq is the unique normalized abelian differential of third kind, i.e. holomorphic on 
£ \ {P, Q} with simple poles at P and Q with residues respectively +1 and —1 and vanishing 
A-periods, while $° = (jPjUJi is a Lie algebra valued holomorphic differential and uj is a basis 
of normalized holomorphic differentials on S. 

As explained in [62] , the Poisson brackets are induced by the Lie algebra of the complexified 
gauge group. At every puncture the residues are expanded as = v^t a ', where t a is a basis 
of the Lie algebra and the Poisson brackets are {z/£, v\ }pb = i5kif abc Vk- 

The Hamiltonians are the Chern polynomials of the Higgs field, namely the coefficients of 
the expansion of the spectral curve (12.11) as 

det ($ - x ■ 1) = 0. 

The quantization of the integrable system, as proposed in [62], is induced by the quantization 
of the Poisson brackets above. 

The case of our interest is a projection to the Cartan degrees of freedom of the general 
integrable system specified to SL(2, C). 

In the sphere case, there are no holomorphic differentials and the Higgs field reads (see also 

TO 

v k dz 



Z — Wk 27TZ 



where v k is the only Cartan element. The corresponding relevant Hamiltonians are generated 
by 

Tr$ 2 = y^ 31 s9 + -^-) (2.4) 
-V V U - w k Y z-Wk \2m J 
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where 



Jl = Tr vl H< 0) = 2 V Tr v,v k (2.5) 



According to the general discussion above, the quantization of these operators is provided by 
replacing v k at each puncture with the corresponding spin operators. 
Analogously, in the torus case, the Higgs field is 

s / (z — w k ) „ . \ dz ,„ „. 

$= 5>*T7 i + 2 ™p)7T- 2.6 

l ^ - w k ) J 2m 



from which it follows that 



Tr* a = Wp(s-u; fc )J2 
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where ' 



H^, ; = 2 > Tr v k v x — + Am Tr u k p, 

Vi{w k - Wi) 



l^k 

k k,l;k^l iV K U 

See Appendix A for the definition of the theta functions. In ( 12. 71) . P is the Weierstrass P- 
function. To obtain ( 12. 7p we used the identity ( 1A.8[) relating the Weierstrass P-function and 
its primitive ('(z) = —V(z). 

In what follows, we will see the appearing of the above Hitchin Hamiltonians via the gen- 
eralized beta-deformed matrix model. 



3 Matrix model: genus zero 

The AGT relation associated with a sphere is the one between the Nekrasov partition function 
of = 2 superconformal SU(2) n ~ 3 linear quiver gauge theory and the n-point conformal 
block on the sphere. Both of them are specified by the marking of Co,n We can obtain the 
beta-deformed matrix model starting from the Dotsenko-Fateev integral representation of the 
conformal block [63] as follows. (See [371 H7J). In terms of the free field <j)(z) (whose OPE is 

t Notice that with respect to (4.10) in [62], in H^ 1 ' we find also a term proportional to r\\ , which will reveal 
to be crucial in comparing with the quantum Seiberg-Witten curve. This term was invisible to the authors of 
[6"2"| . being the absolute normalization of the differential which they admittedly do not check. 
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4>(z)4>(u) ~ — | log(z — 00)), the n-point conformal block is described by inserting the screening 
operators 

Z C - = Uj d\j : e 2b ^ -)j N liv mk (w k )\ , (3.1) 

\ k=0 I f re e on Co 

where the vertex operator V mk (w k ) is given by : e 2rrik ^ Wk ^ :. The momentum conservation 
condition relates the external momenta and the number of integrals as ^fc=o m & = bN — Q. 
By evaluating the OPEs, it is easy to obtain 

Z c ^ = C{m k , w k )Z c ^ = e FC °' n/9 % (3.2) 

where Z c °- n is the beta-deformation of one matrix model 

r N 

Z c ^ = / Y[d\j Y[(\j - Xj)- 2b2 e^ w ^ = e^/rf (3.3) 

"* 1=1 KJ 

and 

n-2 „ 
W(z) = ^22m k hg(z - w k ), C(m k ,w k ) = (w k - w e ) . (3.4) 

k=0 k<£<n-2 

We have introduced the parameter q s by rescaling m k — > — . We relate m k with the mass 
parameters of the gauge theory. Also, we have chosen three insertion points as wq — 0, W\ — 1 
and w n ^i = 00. The remaining parameters are identified with the gauge theory coupling 
constants = e 2ntTi (z = 1, . . . , n — 3) as follows: 

w 2 = qi, w 3 = q 1 q 2 , . .., w n _ 2 = <?i<?2 ■ ■ ■ q n -3- (3-5) 

While the dependence on m n _i disappeared in the potential, this is recovered by the momentum 
conservation condition 

n-2 

y^ j m k + m n ^ x = bg s N - g s Q. (3.6) 

We will refer to F m as free energy. 

The identification of the parameter b with the Nekrasov's deformation parameters is given 

by 

ei = bg„ e 2 = y. (3.7) 

Note that, in the case of b — % (c = 1), this reduces to the usual hermitian matrix model and 
this case corresponds to the self-dual background e\ = —e 2 . 
Here we define the resolvent of the matrix model as 

R(z u ...,z„) = (bg s ) k J2 - E t\~ ■ (3 - 8) 

h h k 

For k — 1, this reduces to the usual resolvent. 
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3.1 Wave-function and conformal block 



In the following, we mainly concentrate on the limit where e 2 — > with e\ and the other 
parameters keeping fixed. In other words, the limit is b — > 00 and g s — > with bg s and N 
keeping finite. This is the limit by Nekrasov and Shatashvili [TT]. In [TQl [29], it was shown 
that the the conformal blocks on a sphere with the additional insertion of the degenerate fields 
V_l(z) = capture the quantization of the integrable systems. 

In this limit, the beta-deformed partition function can be written as J Ylf=i d\i exp(— j-W) 
where 

W = Y1 W(Xj) + 2 ei lo s(A/ " A ^)- ( 3 - 9 ) 

/ KJ 

Thus the leading order part of the free energy can be obtained from the value of the critical 
points which solve the equations of motion: 

W'(Xj) + 2e a V —Up = 0. (3.10) 

We note that these two terms are of the same order in the limit because N and t\ are kept 
finite. 

In this section, we will show that under the identification the beta-deformed matrix model 
Z c °' n with the n-point conformal block, the integral representation of the degenerate conformal 
block can be written in terms of the resolvent of the original matrix model (13.8f) . in the €2 — ^ 
limit. More explicitly, we will show 

^§^(z u z 2 , . . . , z n ) = f[ Wife), %( Zi ) = exp (J- j 1 x{z')dz^j , (3.11) 

in the e 2 — > limit, where Z d °'g +n is the matrix model partition function corresponding to 
the 2n-point conformal block with n degenerate fields inserted at z = Zj. This property of 
"separation of variables" agrees with the corresponding result of the Virasoro conformal block as 
in [32], [10] - The differential x(z)dz is identified with the "quantized" Seiberg-Witten differential 
which is given, in terms of matrix model language, by 

W'(z) 

x{z) = (R(z)) - -^fl, (3.12) 

where (...) = J U d\ [I(A/ - A j)~ 2b2 e b ^ w/9s .... This relation f l3TTD is a simple general- 
ization of the one obtained in [20] for the single degenerate field insertion. 
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First of all, we consider a more generic expression: the (n + £)-point conformal block where 
degenerate fields are inserted 

N n-l 



^r'=(n^w(/^ (A> ) n^w 

\ % — 1 k — t 

TT - 1 TT _ 2m k m e ______ 

= [[(z i -z j ) 3i* (wk-w/,) [(Zi-Wk) 



i<j 0<k<e<n-2 i=l fc=0 

AT n-2 

X 



/_ » / b — ■ 

j=i /<; / fc=o t=i 

where the potential W(z) is the same as (13.41) . We have taken w n _i to infinity and omitted the 
factor including this, as we have done above. The momentum conservation is however modified 
by the degenerate field insertion as 

J2 m k + ^ = bg s N-g s Q. (3.14) 

fc=0 

By dividing by Z Co ' n and taking a log, we obtain 

■<* w = _><* - - D isr + log (n(* - A '>) ■ < 3 - 15 > 

i<j i ys \ i,I I 

Notice that the expectation value is defined as above, but with the modified momentum con- 
servation (I3.14p . By defining e L = YLn( z i ~ we notice that 

/Zi J I 

-Hr' ( 3 - 16 ) 
z i -X 1 

where we have ignored irrelevant terms due to the end points of the integrations. Then, we 
use that the expectation value of e L can be written as log (e L ) = J2T=i M (L k ) conn (2D]' where 
(. . -) conn means the connected part of the correlator, (L 2 ) conn = (L 2 ) — (L) 2 , etc. Thus, (13. 15ft 
can be expressed, by using A3. 16j) . as 

-w-?— + £hUX7 — JL; (317) 

We will now consider the limit where 62 — > 0. In this limit, the terms with k > 1 are 
subleading contributions compared with the k = 1 terms since the connected part of the 
expectation value can be ignored in this limit. Thus, we obtain 



where we have used (I3.12p . Thus, by setting I = n, we have obtained (13. lip . This indicates that 
the properties of the conformal block with degenerate field insertions are build in the resolvent 
of the matrix model in the e 2 — > limit. 



3.2 Loop equations 

The argument in the previous section shows that the relation with the integrable system can 
be seen by analyzing the resolvent, in particular, the loop equations. Thus, we derive it here 
with finite /3. First of all, we keep the potential arbitrary and obtain 







iY 



i=i 



K 



d 

dX 



K 



-h R ^ z ))- b -±i(R( z y) 



zoo 



z- X 



KJ 



9s 9s 9t 



W'(z)(R(z))- 



9l 



where R' is the derivative of the resolvent and we have defined 



(3.19) 



(3.20) 



The expectation value is defined as the matrix model average. By multiplying (I3.19P by g 2 s , we 
obtain 







■(R(z, z)) - ( Cl + e 2 )(R(z)>) + (R(z))W'(z) - f(z). 



(3.21) 



In the case of the hermitian matrix model b = i, the second term vanishes and the equation 
reduces to the well-known one. 

Let us then analyze f(z) by specifying the potential to the Penner type one (13 .4p . In this 



case, 



/(*) 



n-2 



^— ^ Z — Wk 

k=0 K 



(3.22) 



where for k > 2 



Cfc = -bg s 



2mi 



9l 



aiogZ c °." dF n 



dwi 



(3.23) 



While we cannot write Cq and C\ as above because we have chosen wq = and Wi = 1, we can 
see that they are written in terms of Ck with k > 2. First of all, due to the equations of motion: 



W(Xi)) = 0, the sum of Ck is constrained to vanish ^22=0 Ck = ®- m or der to find another 



S 



constraint, we consider the asymptotic at large z of the loop equation. The asymptotic of the 
resolvent is (R(z)} ~ , so that the leading terms at large z in the loop equations satisfy 

n-2 n-2 

- (bg s N) 2 + (ei + e 2 )bg s N + bg s N ^ 2m k - ^ u> fc c fc = 0. (3.24) 

fe=0 k=0 

The leading term of order 1/z in /(z) vanishes via the first constraint. Thus, we obtain 

n-2 /n-2 \ /n-2 \ 

^2 Wk ° k = ( mk + m »- 1 + ifr + ^ ) ( S mfe ~ mn - 1 + ) = A/2 ' ( 3 - 25 ) 

k=0 \k=0 J \k=0 J 



where we have used the momentum conservations (I3.14p with £ = n. Therefore, c and C\ can 
be written in terms of Ck (I3.23p . These constraints are related to the Virasoro constraints [64J. 

€2 — > limit 

As above, in the e 2 — > limit, the connected part of (13. 8p can be ignored in this limit: 
(R(z, z)) — > (R(z)) 2 . Taking this into account, the loop equation (13.211) becomes 

= ~{R{z)) 2 - e 1 {R(z)') + (R(z))W'(z) - /», (3.26) 

where R and / are R\ t2 ^o and f\ t2 ^o respectively. In the following, we will omit the tildes of 
R and /. Then, in terms of x = (R(z)) — w ^ ; the equation becomes 

= -x 2 -e x x' + U(z), (3.27) 

where 

U(z) = —^-jw"(z)-f(z). (3.28) 

This equation is similar to the one obtained in [57]. It is easy to see that this can be written 
as the Schrodinger-type equation: 

= -e 2 1 ^(z) + U(z)^(z), (3.29) 

where ^(2) is defined in (13.111) . 

The above argument is applicable for an arbitrary potential W(z). Here we return to the 
Penner-type one (13 .4p and see the relation with the Gaudin Hamiltoniaril. (I3.28P becomes in 
this case 

= g m fc (m fc + ei ) y, Jh_ _ g 

to ( z ~ W *y k z ~ w * to Z ~ W « 



$We are using the modified momentum conservation to apply this to the argument in the previous section. 
However, the difference will disappear in the £2 — > limit. 

§A relation with Gaudin system at finite £2 has been noticed also in [57] 
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where 

U>fc — We 

U(z) is the vacuum expectation value of (12. 4p . In particular, notice that the residue of the 
quadratic pole in ( I3.30|) corresponds to the eigenvalue of J 2 , quantized in e\ units and that 
H k — Cfc are the vacuum energies of the quantum Hamiltonians (12. 5p . 

Let us rewrite c k in terms of the gauge theory variables. Since the moduli of the sphere 
are related with the gauge coupling constants of the linear quiver gauge theory as in (I3.5p . the 
derivatives with respect to w k can be written as 2iriw k j^- = Q ® — where the second term 
vanished when k = n — 2. Therefore, for k — 2, . . . , n — 2, by using (13. 23j) . we obtain 

1 ( dF Co ' n dF Co - n \ 1 

c k = 7T~ — "a 5 = 7T~ — ~ Uk ) ■ ( 3 - 32 ) 

2mw k y dr fe _i dr k J 2mw k 

where u k are closely related with the gauge theory variables (tr^ 2 ,), <p k being the vector multiplet 
scalar of the k-th gauge group. Indeed, supposing that the free energy F °' n is identified with 
the prepotential (with ei and e 2 ) of the gauge theory (this has indeed been checked for n = 4 
in [551 EH HH1 06] in some orders in the moduli), we can use the e-deformed version [651 EHl [9] 
of the Matone relation [671 EHl EH] to relate u k with (tr^ 2 .). Note that there is still difference 
between F c °' n and F c °' n , we will explicitly consider this in an example below. Instead, for cq 
and Ci, we can use the two constraints derived above and obtain 

n— 2 j n—2 - 

Co = - l)c fc - M 2 = — V ^^(n fc _x - u fe ) - M 2 , 

z — ' 27TZ *— ' W k 

k=2 k=2 K 

n-2 ^ n-2 

Cl = - y WkCk + M 2 = -— YWi - u fc ) + M 2 . (3.33) 

fc=2 fc=2 

We have obtained the differential equations for \I/ = J^f^ ^i{zj) each of which is (I3.29P satisfied 
by ^i(zi). These are similar to the differential equations satisfied by the (n + n) -point Virasoro 
conformal block where n vertex operators are chosen to be degenerate V(z) = e - ^, as in 



|10j . As an example, we will explicitly see in the subsequent section the differential equation 
for n = 4 is the same as that obtained from the Virasoro conformal block. 

3.3 Sphere with four punctures 

We now consider the case corresponding to a sphere with four puncture where the matrix model 
potential is given by 

2 

W(z) = ^2 2m fc ^°&( z ~ w fc)> = 0, Wi = 1, w 2 = q. (3.34) 

fc=0 
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This corresponds to SU(2) gauge theory with four flavors and the Gaudin model on a sphere 
with four punctures where the number of commuting Hamiltonian is just one. Thus we simply 
consider the wave-function ^(z) = e 7 ^ ^ xdz and the differential equation satisfied by it. 
In this case, U (z) can be evaluated in terms of c 2 as 

TT< y \ _ V" m fc + 6imfc , m l ~ E?=o m i + £ i( m 3 ~ ELo m fc) 
[ ) ~^ (z-wuY z(z-l) 



k=0 



-q(q - l)c 2 + 2{qm x m 2 + (q - l)m 2 m ) 
2(2 - l)(z - g) 

Let us relate this with the one obtained from the Virasoro conformal block. We consider the 
last line of U(z). Let us recall the definition of the free energy (13. 2p and (13. 3p . Since the 
difference of them is expressed by C = q 3 (1 — g) »s in this case, the free energies are 
related by 

F c ,4 = pC 0A _ 2mi m 2 l g(l - q) - 2m 2 m log q. (3.36) 
Therefore, its derivative is 

q(l - q)— = q(l - q)— h 2qm 1 m 2 - 2(1 - q)m 2 m . (3.37) 

oq oq 

We notice that the right hand side is the numerator of the last line of U(z). 
Here let us redefine the mass parameters as 

m = m + Y> ^3 = ^3-^-- (3.38) 
In this notation, U can be written as 



7,2 



m, 



U(z) = + 



- -j mi(mi + ei) m 2 (m 2 + e 1 ) 



z 2 (z — l) 2 (z — q) 2 

fh\ + ttlq + mi (mi + ei) + m 2 {m 2 + ei) g(l — g) 



2(2; — 1) z(z — l)(z — q) dq 



(3.39) 



Note here that the first four terms are exactly the potential which considered in [29] V(z). 
Also, the last term might correspond to the "eigenvalue" in [29]. The Schrodinger-like equation 
becomes 

Note that this differential equation has also been derived in [2U] from the free field expression, 
the first line of fl3TT3"j) . 
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4 Matrix model: genus one 



In this section, we consider the matrix model corresponding to the conformal block on a torus 
with punctures [37J EH EEJ EE]. We will derive the loop equations of the matrix model and 
relate it with the differential equations of the corresponding Hitchin system. 

We consider the n-point conformal block on a torus whose integral description is 

-2b 2 n 26mt 



x 



n 



r](r) 

#l(Wfc - Wj) 

r/(r) 3 



#i(A/ - w k ) 



I k=l 

4,TT P Y,k m k w k 



e 9s ^! 1 



= C(w k ,m k ,p)Z Cl ' n , 
where C(w k ,m k ,p) is the A independent coefficient 
C(w k ,m k ,p) =7] 3s s 3 



(4.1; 



ip 

k<i 



i{w k -Wi) b b e 



2m k m e ^vY.k m k w t 
— 1 — 



(4.2) 



and thus, 



Z&,n = e F Cl > n /g 2 a 



/n^n^(^-^)" 262e - E/W(Aj) ( 4 - 3 ) 

J I KJ 



with the potential 



W(z) = ^ 2m k log (z — w k ) + Anpz. 



(4.4) 



k=l 



While we do not know the representation of this in terms of a matrix, we refer to this as 
(generalized) matrix model. As the matrix model in the case of the sphere, this expression 
has been obtained [Ml \5U\ from the free field expression, similar to f l3.ll) but on the torus, 
following from the Liouville correlator by the method in [71]. The parameters must satisfy the 
momentum conservation 



m k = bg s N. 



(4.5) 



fc=i 



See Appendix A for the definition of the elliptic theta functions. 

This matrix model is related with the M = 2 elliptic SU(2) quiver gauge theory [61] which 
is obtained from two M5-branes on Ci.„ with specifying its marking. The gauge theory coupling 
constants q\ = e 2ntTi (i = 1, . . . , n) are identified with the moduli of the torus as 

n— 1 n 

e ^ k = Y[q u q = e ^ = Y[ qi . (4.6) 



i=k 



i=l 
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The parameters m k are directly identified with the mass parameters of the bifundamentals. 
The remaining parameters of the matrix model, the filling fractions and momentum p in the 
potential determine the Coulomb moduli of the gauge theory. 



4.1 Wave- function and conformal block 



As in the previous section, we can show the relation between the resolvent of this matrix model 
and the conformal block with the degenerate fields $1^ = e~^( z ). The integral representation 
of the latter is 



J I KJ I 



2bm 1 4wbp 
Ss g 9s 



(4.7) 



n n 



i=l 



r)(r) 3 



n 



_ 2?T£, 



n 

Kj L 



77(7-) 



1 

" 2p 



where Zi are the insertion points of the degenerate fields and the momentum conservation is 
slightly deformed from the original one to 



fe=i 



J2 mk + ^ = b9sN - 



(4.{ 



7 C l,n+l 



As before, we consider log p? ln . The main object we want to know is log \ i $i( z i ~ Aj) », 
where the expectation value is defined in the same way as in the sphere case. Note that we 
are using the momentum conservation (14.81) . Then, if we define e L = Yli i$i( z i — A/), we 
can rewrite it as L = J2i i By using the same argument as the sphere case: 

log (e L ) = Efc n ( Lk ) CO nn' we tnerefore obtain 



log 



J deg 



^W{ Zi ) 31(1+1) , 1 ^ 

% yA i<j 



k=l 



(4.9) 



where we have used the momentum conservation in the second term. 

We introduce the deformation parameters and take the limit where — > while e\ keeping 
fixed. In this limit, the path integral is dominated by the solutions of the equations of motion 

^i(Aj-Aj) 



W'(Xj) - 2b 9s J 



3& 



0. 



(4.10) 
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as we have seen in the previous section, and the connected part with k > 2 in (14. 9 p is negligible. 
Thus we obtain 

Z Cl ' n+e i r z i / W(?'\\ 

log^^fE/ ((^))-^)^-iogn^)- (4.H) 

^ i ^ ' i 

where R is an analog of the resolvent of the sphere case 

r( Zu ..., Zk) = dg.)* t'f'~M ■ ■ ■ e < 4 - i2 > 

^ i>i(zi - \ h ) ^ Vi{Zk - A/J 

Thus, we have related Z d ^" +t with the resolvent or ^ at least in the 62 — >• limit. 

One simple consequence which follows from the above formulas is about the monodromies 
of or Z de 'g 1+1 . (For simplicity, we consider the £ — 1 case. The generalization to £ > 2 might 
be straightforward.) Along the A cycle, Z de ^ x behaves as 

+ «) = i-lf-^e-^^iz) = (-l)^e-^% C e 7 +1 (,), (4.13) 

where we have used the momentum conservation. Similarly, we can evaluate the S-cycle mon- 
odromy 

Z deq +1 ( Z + 7rT ) -p 3 /rr« / t , 2iV,A,\ ( %Z WlT ^Tp 



e-^trj-*? ( ~ A 7 )e 2l ^M exp - — - -^J (4.14) 



On the other hand, let us consider the shift of p — > p — in Z d ^ +1 , which gives the additional 
factor exp (— |f + 2i J2i ^i) m t ne integrals. Therefore, we obtain 

We note that these are indeed the same monodromies as those of the conformal block with the 
degenerate field [301 123- 

We can further proceed to derive the special geometry (ei-deformed Seiberg-Witten) relation 

Cl, n 

deg 



for the resolvent. In the e 2 — > limit, (14. lip shows that Z d ^" +1 can be expanded as 



Z^; + \z) = exp + I J Z X ( z ') dz > + (e 2 )^ , (4.16) 

where the first term is ^(ei) = lim^^o F Cl ' n . The above y4 and S cycle monodromies indicate 
that for the second term 

x(z)dz = -2tt 2 p, / x(z)dz = -2tx 2 tp - 1 dJ7 ( 6l \ (4.17) 
Jb 21X1 d P 

This corresponds to the Seiberg-Witten relation in the presence of the e\ dependence [121 13] • 

For other independent cycles which correspond to the legs of the pants decomposition of the 

torus, it is natural to expect that the similar relation is satisfied. 
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4.2 Loop equations 

Let us derive the loop equations of the matrix model for the torus. As discussed above, in order 
to relate with the conformal block, we will use the momentum conservation ( 14. 8ft . From the 
Schwinger-Dyson equation for an arbitrary transformation 5Xk = , we derive 

*-* (? (M^f) - * (? mt) + ^« (? 

where we have multiplied the both sides by g 2 and defined 



t{z) = bg s ^[l^ iW'M - W(z))J . (4.19) 
We then use the formula (1A.9|) to calculate the last term and, after some algebra, we obtain 
= - (R(z, z)) - [b + ^ s <#(*)) + W(z) (R(z)) + 6 2 ^ 2 iV ^ H^f ^ + *(*) 

(E Sa'Ia') ) + 3b2g2 ^ N ~ ^ ^ 

This equation is valid for an arbitrary potential. 

From now on, let us consider the potential corresponding to the toric conformal block (14. 4p . 
In this case, t(z) can be evaluated by using (1A.9|) again as 

±f \ 01 &'i(z ~ Wk) /v^ #i(A/ - Wk)\ h /K( x i- W k) 

t{z) = 2b9 ° E mk M^) \E ^T^y ) - <*• E ^ \ ^(a7^) 

-^iV £ m fc f r (Z " W \ } - ^ E ™* (E f!""M \ - E ( 4 - 21 ) 

By substituting this into (14. 201) and using the momentum conservation, we obtain 

= - (R{z, z)) - (b + ij <? s (#(2)) + - 3fo^(iV + 1)^ ^ m k 

+ 26 *' E mfc X7T— TT E ^7Y73^ - 6 ^ E m * 



fe=i 



' $i{z - w k ) y-f ftiiXi - w h ) tf^z-Wk) 



-^E^(E^7^/ +6 ^\g^A^ 

*-A 7 ) 
-At) 



git ^€{z-Xi)\ U . 
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where the last two terms come from the deformation of the momentum conservation. Let us 
note that the r-derivative of the partition function Z Cl ' n 

Z^n ding ^2. m * \2_. MXl - Wk ) I + b9s i?i(Aj — Aj) y ' (4 - 2dj 

where we have used — 8^-i9i(a) = This is the third line in ( I4.22p . Similarly, the 

derivatives with respect to w k produce the first term in the second line of (14.221) : 

- 2bg s m k ( V ^ . (4.24) 



We will fix one of the moduli of the torus as w n = below. In this case, the derivative with 
respect to w n is understood as the right hand side of (I4.24p . Putting all these together, we 
finally obtain 

= - (R(z, z)) - (b + ij g s (R'(z)) + W\z) (R(z)) - 3bg s (N + 1)^ £ m k 

o^VUz -w k )d\nZ c ^ , ArV ^ $Uz-w k ) A 2 d\nZ c ^ 
£^#1(2 -■">*,) <9w fc ^ V^z-Wk) d\nq 



e2 — > limit 

We consider the e<i — > limit. We can use under this limit the equations of motion ( I4.10p . 
Furthermore, as discussed in the previous section, (R(z,z)) ~ (R(z)) 2 . Also, the last two 
terms in ( I4.25P disappears in this limit. Naively, the terms with the derivative of \nZ Cl - n also 
disappear due to the factor g 2 s {= e^). However, they does not because lnZ Cl - n behaves as 
F Cl >™ i 1 g 2 s . Thus, the loop equation (14.251) becomes 

= - (R(z)) 2 I - ei (R'(z)) I + W\z) (R(z)) - 3 Vl (V ™t)(Y] m k + ei) 



I k 



V m^)> m k -±) 4-}^) ^-^ +4^ 4.26 



A 



$i(z-w fc ) ^^l^-^fe) dw fc 



£2^0 ding I e 2 ^o 



We will omit in what follows. After some algebra (by introducing x = {R(z)) — W/2 and 
by using the formula of the theta function), we obtain 

= -x 2 - eix' + ^2 m k(m k + e 1 )V(z - w k ) 

k 

^ Vi(z - w k ) y ow k J omq 
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where V(z) is the Weierstrass elliptic function (1A.3j) and 



H k = 4:7ipm k + 2 



m k m t - 



H = 4vrV - m V m k [m k + e x ) + J V m,m j ~ ^ . (4.28) 

Therefore, we obtain that H k — d ^' n and H + 4 are the vacuum energies of the quantum 
Hitchin Hamiltonians Hjj. and Hq 1 ^ in ( 12.7ft . This shows that this matrix model captures the 
quantization of the Hitchin system associated with the torus. 

We can write this equation in the form of the differential equation satisfied by the wave- 
function ty(z) = e~ f x ^ z ^ dz : 

<9z 2 ^ ^ - W k ) J 

By considering n^Li ^(^i)^ this satisfies the above differential equations for each Zj. This is 
related with the KZB equation [72j [73j [74] as discussed in [28] (see also [75]). 

It is easy to translate the terms in the right hand side in (14.291) to the gauge theory variables. 
By the identification of the moduli, the derivatives with respect to w k are written as = 
■Jr jj-^ — (for k — 1, . . . , n — 1), where tq = r n , and also t^t— = tt^jt-- Therefore, we obtain 

OTfc OTfc-i v 11/1 v in ding 2m ot„ ' 



dw k 



u k - u k -i, (4.30) 



for k = 1, . . . ,n — 1, where we have defined u k = 9F d ^ k n arid u$ = u n . For k = n, we calculate 
^ = - > , h 4vrp 2_^m k = -u n ^ +u n + Airp ^ m fc (4.31) 

n fc=l k k=l k=l 

where we have used the equations of motion and the momentum conservation. (We also ignored 
the term depending on 62-) Also, for the derivative with respect to q we obtain 

- 4—— = : u n . (4.32) 

a in q m 

As in the sphere case, these u k are related with (tr<^). Note however that there could be a 
difference of them since u k here are the derivatives of F, which will be seen below. 
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For completeness, let us rewrite the above equation in terms of the original partition function 
Z Cl ' n . The difference between Z Cl - n and Z Cl - n is given by C (14.21) . which gives rise to 



=4^j h 3 ) m k (m k + ei)?7i - ) m k m e — -, 

alng c>lng |— ^ Vi(w k — Wi) 

— a = q 2 } m k m e — - Anpm k , (k = 1, . . . , n - 1 

dw fc ^ wi(ty fc -u;/) 



n fc=l 16 f(^ra) V y 

where we have used (14.311) . Thus, we obtain from (I4.27P 
= -x 2 - exx' + ^2 m k(m k + e x )V{z - w k ) + 2 ^ m fc (m fc + e^i (4.34) 

fc k 

_ ^ y l( .- mt )a^ + + y W\ + 4lr y + 4^. 

~0i(*-iUk) <9u> fc #1(2) 1^ <9w fc ^ $x(we) I d\nq 

4.3 One- punctured torus 

In the n — 1 case, we can see the relation with the elliptic Calogero-Moser model. We will take 
wi = 0. The potential is 

W{z) = 2mx log 1 (z) + 4vrpz. (4.35) 
In this case, it is easy to calculate the loop equation (I4.34j) : 

= -x{z) 2 - exx'(z) + mx(mx + ex)V(z) - 4u(ei), (4.36) 

where 

, , 22 ndlnZ 01 ' 1 mx(mx + ex)r}x 

U ^ = ~*P -9°^q 2 

d 

= -*V + 771 — [J* 1 ' 1 ~ 2mx(mx + e t ) In 77) , (4.37) 
amg 7 

where we have used that rjx = 4§{^ (See Appendix A) and defined the free energy as J rCia = 
lim^-^o F x > x . Note that the free energy is the one evaluated in the 62 — > limit. Eq. (14.371) is 
the ei-deformed version of the relation between (tr0 2 ) and the Coulomb modulus u [66] which 
coincides with the one found in |24j . 

By introducing the "wave- function" \1> = e e i , we finally obtain 

d 2 



^ = 4u^. (4.38) 
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The left hand side is the Calogero-Moser Hamiltonian and u in the right hand side can be 
considered as the eigenvalue of the Hamiltonian. Note that similar equations have been derived 
from the Virasoro conformal block [29] and affine conformal block [28]. Indeed, by assuming 
the equivalence of the partition function of the matrix model Z 1,1 with the one-point conformal 
block on a torus, ( 14. 38p becomes the exactly same equation as the one obtained from the 
conformal block with the degenerate field. (See Section 3.1.2 in [29]. The identification of the 
parameter is a = inp.) 

We also emphasize that the differential xdz in the wave-function satisfies the special ge- 
ometry relation ( 14.1 7ft . This is equivalent to the proposal in [12] stating that the ei-deformed 
prepotential can be obtained from the ei-deformed special geometry relation for the M = 2* 
theory, by using the same argument as in [29] . 

4.3.1 Large N limit and prepotential 

Before going to next, let us consider the loop equation in the large N limit which can be obtain 
by taking e± — > further in (14.361) : 

x 2 = m 2 P(z) -4«(ei = 0), (4.39) 

which is the Seiberg-Witten curve of the Af = 2* gauge theory [3]. Indeed, the parameter u 
can be written as 

u(ei = 0) = -Try + TT, — ( tf 1 ' 1 ~ 2m l l n v) ■ (4-40) 

amg V / 

where Jo' 1 is the leading contribution of the full free energy in the limit where ei^ — > 0. The 
first term corresponds the classical contribution to the prepotential. The last term denotes the 
shift of the Coulomb moduli parameter from the value of the physical expectation value (Tr </> 2 ) 

m Esi El]. 

Indeed, we can be more precise. Under the identification inp with the vev of the vector 
multiplet scalar a, it is easy to show from f 14 . 1 3 j) that 

i r 

a = inp = - — : / xdz. (4-41) 

This and the fact that the form of x here is the same as the Seiberg-Witten differential of 
the Af = 2* gauge theory where J^ 1 ' 1 is changed to the prepotential (see [771 122]) lead to the 
conclusion that the free energy in the large TV limit of this matrix model is exactly the same as 
the prepotential of the gauge theory (under the identification above). 
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5 Conclusions 



In this paper we proposed that the /3-deformation of matrix models provides, in a suitable limit, 
the quantization of the associated integrable system. In particular we have shown that the loop 
equations for the /3-deformed generalized matrix models [541 [56] reproduce in the Nekrasov- 
Shatashvili limit the Hamiltonians of the quantum Hitchin system associated to the sphere and 
torus with marked points. Moreover, we have shown how to obtain the wave function from 
degenerate field insertions. 

It would be interesting to understand if this procedure could provide a general quantization 
prescription of integrable systems which can be linked to specific matrix models. To this end 
it would be very useful to provide further evidence and examples. For instance, it would be 
interesting to investigate the /3-deformed Chern-Simons matrix model jl9] in this direction. 
Furthermore, the extension of our approach to q-deformed conformal blocks, along the lines of 
[78] , would be worth to be analyzed with the aim of connecting our results with topological 
strings. 

On a more specific side, a natural extension of our analysis concerns Hitchin systems on 
curves of higher genera, the point being a generalization of the identity (]A.8|) . A further 
explorable direction would be the extension to higher rank gauge groups with a multi-matrix 
model approach. 

The problem of understanding the proper quantization of the Seiberg-Witten geometry has 
been explored recently also from a different view point consisting in a saddle point analysis of 
the instanton partition in the Nekrasov-Shatashvili limit [191 [2H [27J . I n Section H] we have seen 
that for the M = 2* theory the eigenvalues of the Hamiltonian match. It would be interesting 
to further explore the relation between the two quantizations. 
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Appendix 



A Elliptic functions 

The elliptic theta function is defined by 

oo 

$i{z\t) = 2q 1/8 smz -<f)(l - 2q n cos2z + q 2n ) (A.l) 

n=l 

which has pseudo periodicity 

■& 1 (z + tt\t) = -'& 1 (z\t), A(z + ttt\t) =e- i{2z+nr) A(z\T). (A.2) 

This function satisfies i9^(z\t) = — 8-q^$i(z\t) where $[(z\t) = ^i?i(^|r). 

The Weierstrass elliptic function V is double periodic with periods it and tit and is expressed 

as 

V(z) = -C(z) t CW = ^ + ^, ^ = -6 ^|r)U =0 ' (A ' 3) 
where i?i(z|r) is elliptic theta function. The Weierstrass function satisfies 

P(^)' = 4P( 2 ) 3 - s 2 7^) - S3, (A.4) 

where 



#2 : 



4 / v-^ v? q n 

- 1 + 240 > — 

3 I ^ 1 - q n 



q 

n=l 



oo 



S3 = i(l-™E^J, (A-5) 
We also define gi whose expansion is 



q' 

n=l 



3C 



n=l 



which is related with rji as g\ = —2rji: 



i?i(6-a) V^i(a) #1(0)^(6) 2 \^(a) ^(6) 0i(6-a) 



(A.7) 



where r\ = g 1 / 24 I3^ =1 (l — g") is the Dedekind eta function. 
The Weierstrass and zeta functions satisfies the identity 

V(a + b)+ V(a) + V(b) = (C(a + b) - C(a) - C(&)) 2 - (A.8) 

This leads to the following formula of the theta function 

^{b-a) (d'M ^{b)\ ^(a)^(b) 1 f^(a) , , 0?(6-a)^ 
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